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Abstract 



We consider in this work a model for aggregation, where the coalescing particles ini- 
tially have a certain number of potential links (called arms) which are used to perform 
coagulations. There are two types of arms, male and female, and two particles may co- 
agulate only if one has an available male arm, and the other has an available female 
arm. After a coagulation, the used arms are no longer available. We are interested in the 
concentrations of the different types of particles, which are governed by a modification of 
Smoluchowski's coagulation equation — that is, an infinite system of nonlinear differential 
equations. Using generating functions and solving a nonlinear PDE, we show that, up to 
some critical time, there is a unique solution to this equation. The Lagrange Inversion 
Formula allows in some cases to obtain explicit solutions, and to relate our model to two 
recent models for limited aggregation. We also show that, whenever the critical time is 
infinite, the concentrations converge to a state where all arms have disappeared, and the 
£ — . | distribution of the masses is related to the law of the size of some two-type Galton- Watson 

tree. Finally, we consider a microscopic model for coagulation: we construct a sequence 
of Marcus-Lushnikov processes, and show that it converges, before the critical time, to 
the solution of our modified Smoluchowski's equation. 
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1 Introduction 

In this work, we study a model for coagulation of particles, generalizing the original model 
of Smoluchowski [26], and a recent model of Bertoin [2]. We consider particles which are 
initially given a certain number of male and female arms. These arms are used to perform the 
coagulations: two particles coagulate when a male arm of one and a female arm of another bind. 
This can be used to model the formation of polymers. For instance, consider male particles 
(which have only male arms), and female particles. Then a coagulation between a male and 
a female particle can be thought of as an ionic bond between a cation and a anion. This 
kind of models has also been investigated in the physical literature. For instance, in [13], (14] . 
the authors study coalescing monomers with two types, A and B, with bonding only allowed 
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between A and B, hence forming alternating linear polymers. In this work, this corresponds to 
giving to each particle exactly one male arm and one female arm. 

In our model, a particle is characterised by a triple (a, b, m), a G N being its number of male 
arms, b G N its number of female arms, and m G N* its mass. Two particles may coagulate 
when one has an available female arm and the other has an available male arm, and when a 
coagulation occurs, the used arms disappear. Hence, we may only observe the transition 

{(a, b, to), (a', b', to')} — > (a + a' — 1, b + b' — 1, to + to'). 

We will assume that this transition occurs with a rate given by the number of pairs formed of 
a female arm and of a male arm, that is a'b + ah' . We wish to study how the concentration of 
each type of particle evolves when time passes. The precise mathematical formulation is given 
in Section 2. 

This model is a modification of the well-known model of Smoluchowski [26]. Recall that 
Smoluchowski's coagulation equations [26] describe the evolution of the concentrations of parti- 
cles in a medium, where particles are characterised only by their masses. When two particles of 
masses to and to' coagulate, they merge into a single particle of mass to+to'. Such a coagulation 
occurs with rate k(to, to'), where k is some symmetric nonnegative kernel. In Smoluchowski's 
original model, the masses are assumed to be positive integers. The concentration q(to) of par- 
ticles of mass to is governed by the following infinite system of nonlinear differential equations 

m— 1 +oo 

Ct{m!)ct{m — to')k(to', to — to') — q(to) ct(m')K(m, to'), 

m'=l m'=l 

for to G N*. The first term accounts for the creation of particles of mass to by coagulation 
of particles of mass to' and to — to'; the second for disappearance of particles of mass to by 
coagulation with other particles. 

For general kernels k, explicit solutions are not known. However, some have been obtained in 
different cases, notably whenever the kernel is constant [26] , additive [11] or multiplicative [20J . 
In the multiplicative case, solutions are obtained up to a critical time, known as the gelation 
time. This is interpreted as the time when a particle of infinite mass appears. It absorbs some 
of the particles, and the total mass starts to decrease. 

Smoluchowski's equation (and some variations) have been extensively studied, both from 
an analytical (e.g. [5], [8], [15]) and a probabilistic point of view (e.g. [Bj, [16] , [IE], [21], and 
see also the review by Aldous [lj). In general, little is known after the gelation time, and most 
results are obtained before (see however [9], [TO])- The existence and uniqueness of a solution 
before gelation has been obtained only in 1999 by Norris [21], under the assumption that k is 
sublinear, i. e. K(n,m)/(nm) is bounded. 

From a probabilistic point of view, some microscopic models have been studied, beginning 
with Marcus [23] and Lushnikov [2T] . Heuristically, one considers a finite number of particles, 
and each couple of particles with masses to and to' coalesces with rate k(to, to'). After suit- 
able change of time and renormalization, one expects this system to converge to a solution of 
Smoluchowski's equation. This has been shown by Jeon in 1998 [18j (up to extraction of a 
subsequence), provided the rate is strictly sublinear (i.e. K(n,m)/n — > when n — > +oo). In 
particular, there is no gelation in this case. Norris [21] extended his results one year later by 
showing the convergence of the model before gelation, whenever the rate is sublinear. Other 
points of view are also considered; e.g. in [16] . the authors show that coagulating Brownian 
particles form clusters whose size evolves according to Smoluchowski's equation. 




2 



An interesting question is to deal with the case when the coagulations are restricted by some 
device. Typically, one may think of covalent bonds: a given atom can only perform a given 
number of bonds. In this direction, Bertoin [2] studied two models where a particle is charac- 
terised by its number of arms and by its mass, and it uses its arms to perform aggregations. 
The concentrations of each type of particle is governed by a modification of Smoluchowski's 
equation. In [2], he obtains solution up to some time T, and shows that whenever gelation 
does not occur (i.e. T = +00), there is a limit state where all the arms have disappeared: the 
concentrations converge to limiting concentrations which bear a striking resemblance with the 
law of the size of some Galton- Watson tree. This fact is explained in [3] and [I]. It is also worth 
noticing that Bertoin's model can be related to Smoluchowski's for the constant, additive and 
multiplicative kernels. We will also see that our sexed model contains Bertoin's: the oriented 
model corresponds indeed to ours if each particle is given precisely one female arm, and the 
symmetric model corresponds to the sexed one if the particles are given a gender uniformly at 
random. 

This paper is divided in two parts. In the first one (Sections 2 to 5), we shall study the sexed 
Smoluchowski's equation, which is an infinite system of nonlinear differential equations. We first 
(Section 2) introduce the problem, and prove some physically intuitive facts. Then (Section 3), 
we prove our main result: up to some critical time, there exists a unique solution to the system, 
and its moment generating function can be expressed explicity in terms of the initial data. 
The tools used are analoguous to those in [2] , but since we are dealing with a two-dimensional 
problem, several technical issues need to be addressed. The outline of the proof is as follows. 
First, we transform the system into a PDE problem by considering the generating functions of 
the concentrations. This PDE is not quasilinear, but it may however be solved by the method 
of characteristics. This method requires the inversion of a two-dimensional mapping, and this 
can be done precisely up to the critical time. Unfortunately, even for monodisperse initial 
conditions (i.e. there are only particles of mass 1 at time 0), the inversion is not explicit (one 
could use the two- variable Lagrange inversion formula, but in general, the expression it provides 
is too cumbersome). Nonetheless, in some specific cases (Section 4), the Lagrange Inversion 
Formula yields explicit results. In particular, we recover the solutions obtained in [2J. Finally, 
we show (Section 5) that there exist limiting concentrations when t — > +00, and that they are 
related to the distribution of the total progeny of some two-type Galton- Watson process. 

In the second part (Section 6), we study a microscopic model. Given a finite number of 
particles, we let them coagulate and observe the evolution of the concentrations of the different 
types of particles. This is a Marcus-Lushnikov process, and we show that it converges, before 
the critical time, to a process solving Smoluchowski's equation ([1]). As pointed out earlier, this 
kind of convergence had already been proved by Norris ( |23] , see as well [18] ) . The difference 
here is that we consider a model with male or female arms. Moreover, the proof is made much 
easier by the fact that the rate of coagulation is explicit. In particular, we will appeal to 
the PDE obtained in the first part. This discrete model provides a justification to the sexed 
Smoluchowski's equation ([I]). 

Finally, note that our construction can also provide a model for random oriented graphs, 
called the configuration model, since a coagulation can be seen as the creation of an oriented 
edge between two vertices in a graph, whose orientation is given e.g. from the male arm to the 
female arm. Hence, we can consider a large number n of particles and let them coagulate. When 
all the coagulations are performed, we obtain a set of oriented graphs. When n — > +00, we may 
wonder what the distribution of their sizes is, what a typical graph looks like, etc. A heuristic 
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answer, motivated by the works [3], [I], and by the results obtained in this paper (Section 
5), is that a typical graph would be a two-type Galton- Watson tree (with the convention of 
orientation above), provided there are few arms (with the notations of this paper, this means 
T c = +00 and fi is not degenerate). 

2 Setting and results 
2.1 Notations 

Let us first introduce some notations and Smoluchowski's equation, and state our main result. 

• N = {0, 1, 2, . . . } and N* = {1, 2, . . . }. 

• S 1 = N x N x N* is the set of the different types of particles. A generic element of S will 
be denoted by p, and if p = (a, b, m), we will call a p-particle a particle with a male arms, 
b female arms, and mass m. 

• For p = (a, b,m) G S and p' = (a', b', m') G S, we will denote 

p.p' = a'b + ab' 

the rate of coagulation and 

pop = (a + a' — 1, b + b' — 1, m + m') 

the type of the particle resulting from such a coagulation. We say that p' z< p if a' < a + 1, 
b' < b + 1 and ml < m — 1 . When p' -< p, we write 

p\p' = (a + 1 — a', b + 1 — b', m — m!) 

the type of particle such that p' o (p\p') = p. 

• For two functions c, / : S — >• 1R, we will denote, when the series converge absolutely, 

(c,/>:=5>(p)/(p). 

pes 

When using this notation, we will write, with a slight abuse of notation, a for the function 
(a, b, m) 1— > a, 6 for (a, 6, m) — > 6, etc. 

Let us recall our goal. We are interested in a system of coagulating particles with male and 
female arms. We assume that each couple formed of a p-particle and of a p'-particle coagulate 
with rate p.p', to form a p o p'-particle. This means that if we denote Ct{p) the concentration 
of p-particles then (ct(p),p G S) solves the following infinite system of nonlinear differential 
equations 

^ C *(P) = 2 ' ■{P\p') c t{p') c t{p\p') ~ ct(p)^2p-p'c t (p). (1) 

pMp p'es 

The first term accounts for the creation of p-particles by coagulation of p'- and p\p'-particles 
(the factor 1/2 comes from an obvious symmetry). The second accounts for the disappearence 
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of p-particles by coagulation with other particles. Let us once write down this formula explicitly. 
For all (a, b, m) G S, the concentration of (a, b, m)-particles verifies 

, 1 m— 1 o+l 6+1 

l Ct (a, b, m) — E E 5>'( 6 + 1 - & + + 1 - af))x 

rn'=la'=Ob'=0 

Ct{a! , 6', m')c t (a + 1 — a', 6 + 1 — b' , m — m) 
— y^l ^ + a 'b)c t (a, b, m)ct(a', b', ra). 

m'>X a'>0 b'>0 

Let us now define what we call a solution to Smoluchowski's equation. 

Definition 1. We call a family {{ct(p)) p es,t G [0, T)) of differentiable functions a solution of 
Smoluchowski 's equation ( or of system ([I]) ), if 

1. For every t G [0, T), (a + b, \c t \) < +oo, 

2. (a 2 + b 2 , \c t \) < +oo for t in a neighbourhood of 0, 

3. The family (c t (p)) solves the system ([T]) fort G [0,T). 

Remark 1. • We will always assume that at time 0, (a + b + l,Co) < +oo, and that the 
mean number of male arms (a, Co) and the mean number of female arms (b, Cq) are equal. 
Physically, it is then obvious that they will remain equal as time passes. This shall be 
proven later on, in Lemma 

• It is easy to see that if (c t ) te [ 0j T) is a solution to flTJ with initial conditions c , and A > 0, 
then (Xc t /x2)te[o,T) is a solution to flTJ with initial conditions Ac . Hence, it is enough to 
assume that (a, cq) = (b, cq) = 1, what will always be the case from now on. 

2.2 Main result 

Our main result is existence and uniqueness of a solution to flTJ) up to a critical time. In 
all the statements and proofs, we are given nonnegative initial concentrations Co such that 
(1, Co) < +oo, (a, Co) = (b, Co) = 1 and (a 2 + b 2 , cq) < +oo. We can then define the critical time 
T 

Definition 2. Let 

M = (ab, c ) + a/ (a 2 - a, c ) (b 2 - b, c ) 

and 

+oo if M < 1 
ih if M > 1. 

We will also constantly use the generating function of (c ) 

g (x,y,z):= ^ c (a,b,m)x a y b z m . 

(a,b,m)£S 

Since (l,c ) < +oo, g is well-defined on [0, l] 3 . Using the assumption (a + b,c ) = 1, and 
e.g. monotone convergence, we also see that its partial derivatives with respect to x and y are 
well-defined and continuous on [0, l] 3 . For the same reason, they remain in [0,1]. We shall 
prove the following result. 
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Theorem 1. (i) Smoluchowski's equation (Op) with initial conditions Co has a unique solution 
(c 4 ) defined on [0, T c ). 

(ii) For t G [0, T c ), (a 2 + b 2 , c t ) < +oo, and (a 2 + b 2 , c t ) —> +oo w/ien t —>T C . 

(iii) For t G [0, T c ) and 2 G [0, 1], the mapping (p t (-, •, z), given for (x, y) G [0, l] 2 by 

4> t (x, y, z) = ^(i + t)x - t-jj^ixi y> z )i( 1 + t)v ~ f> z ) 

has a right inverse h t = (hi , h^f 1 ) which is well-defined and analytic on (0, l) 2 . Then the 
generating function g t of (q) is given by 

9t(x,y,z) = {H ( t 2 \x 1 y 1 z) + H[ 1 \^y,z))+G t (z). (2) 

where for t > 

fc,f> = l±**W(x, v ,,)-| ; fc»:=l±*ft»(,, B ,,)-J (3) 

and 

• is the antiderivative of with respect to y, vanishing at y = 0, 

~ (2) ~ (2) 

• HI is the antiderivative of h t with respect to x, vanishing at x = 0, 

• Gt(z) is the antiderivative of 



dg 



tf ) (0,0,2),/4 2) (0,0,*),*) (4) 



dz 

with respect to z, vanishing at 0. 
(iv) The total mass (m,ct) is constant on [0,T C ). 

2.3 Preliminary results 

In this section, we give some physically intuitive results, and deduce the "weak" form of the 
equation. Let us start with the following lemma (recall that Ct(p) is meant to model a concen- 
tration) . 

Lemma 1. Any solution to Smoluchowski's equation remains nonnegative, i.e. if ' (ct)te[o,T) ^ s 
a solution to then for all t G [0, T) and p G S, c t {p) > 0. 

Proof. Take some t G [0, T). System (pQ) gives 

^c t (a, b, 1) = -ct(a, b, 1) ^ ^ X^ a& ' + a ' b ) c t( a ' '> b '> m ') := _ 7(£) c t( a , &, !)■ 

m'>la'>0fe'>0 
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Let G(t) = Jq'j(s) ds. Then c t (a,o, 1) = c (a, b, l)e~ G ^\ so it remains nonnegative. Let now 
m > 1, and suppose that the Ct(a,b,m') are nonnegative for a, 6 > and 1 < m' < m. For 
some p = (a, 6, m + 1), we have 

p'<p p'gs 
=f3(t)-c t (ph(t). 

So we may write 

ct(p) = (co(p) + jy^)e G(s) ds) e - G <*>. 

But > since it is a linear combination with nonnegative coefficients of the Ct(a,b,m') 
for a, b > and 1 < w! < m. So c t (a,b,m + 1) is nonnegative, what gives the result by 
induction. □ 

The following two lemmas are straightforward generalizations of Lemma 1 and 2 in [2] . Note 
however that the monotone convergence used in the proofs requires that the coefficients (q) be 
nonnegative. 

Lemma 2. (i) If (c t ) is a solution to Smulochowski's equation then t \— > (l,c t ), t i— > (a, c t ) 
and 1 i — > (6, c t ) are decreasing. 

(ii) A family (c t ) is a solution to ([I]) if and only if it solves 

J t ^f) = \ E P-P' c t{p)ct{p'){f{pop')-f{p)-f{p')) (5) 

for every bounded f : S — > R. 

Remark 2. • The derivative in this lemma has to be understood in the weak sense, i.e. the 
formula actually holds in the integral form. But if f(a,b,m) — > when (a,b,m) — > oo, 
then it is easy to check that the formula holds in the strong sense. 

• Consider in particular, the generating function of c t , g t (x,y,z) = {x a y b z m ,c t ). Then 
g(., ., z) is regular, in the sense of Definition [3] below. 

Definition 3. We say that a function (t,x,y) h- > g t (x,y) defined on [0,T) x (0, l) 2 is regular if 

• 1 1— > y) C 1 and (x, £/) i— > ~j^{ x i V) are C 1 , 

• (x,y) h-> g t (x,y) is C 2 , t h-> —±( x ,y) and 1 1-> -^(x,y) are C 1 . 

Lemma 3. Let (c t ) fre a solution to Smoluchowski's equation, and let 

T r = inf {t > 0, (a 2 + o 2 , q) > r} and = sup r r . 

r>0 

Consider the mean numbers of male and female arms A t = (a, q) and B t = (b, q), and assume 
Ao = B = 1. T/ien 

A = B t = — *— (6) 
1 +t w 

/or a// 1 G [0,T AToo). 
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3 Proof of the theorem 



3.1 Overview of the method 

In this section, we give a sketch of the proof which contains all the important ideas. The 
rigorous proof however requires some care, and it is given in detail afterwards. So, consider a 
solution (ct)te[o,T) to Smoluchowski's equation (CQ), and 

g t (x, y, z) = (x a y b z m , *>=£££ c t (a, b, m)x a y b z m . 

a>0 b>0 m>l 

Using (jHJ) and Lemma El it is easy to see that g t solves the following PDE 

9gt = dg t dg t 1_ (dgt + y ^£t\ 

dt dx dy 1 + 1 \ dx dy J 

Now, we can solve this PDE using the method of characteristics: we want to find a trajectory 
(x(t),y(t)) starting from some (x, y) G [0, l] 2 such that g t (x(t),y(t), z) is easy to compute. So 
let 

(pi(*),P2(*)) = (^«t),y(t),z),^(x(t),y(t),z)\ . 
An easy calculation shows that 

Pl{t) = w { x(t) +P2(t) -—t) + dx^y- [v®+n® ~T+i) ~ IT? (8) 

and a similar formula for p 2 . Now, if we require 

x(t)+p 2 (t)-p£- = y(t)+ Pl (t) //(/) " 



1 + t * w " w 1 + t 
then 

These ODE's are readily solved, with pi(0) = -——{x, y) and WO) = -7— (x, y), and we obtain 

ax dy 

M = ffl (9) 

and 

x(t) =x+(x-p2(0))t ; y(*) = y + (y-Pi(0))*. 
Using the PDE, we now see that 

—gt{x{t),y(t), z) = II y , (10) 

so by integrating 

9t(x(t), y(t), z) = g t ((j)t(x, y, z),z) = g (x, y, z) - 2/' ^O^pfo y ' ( n ) 
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To obtain g t , it only remains to invert <f) t , for, if <f>t(ht) = Id, then 

g t (x, y, z) = g (h t (x, y, z)) - - ^> z^y^fo J/> 2 ))- 

We may now start a rigorous proof, which consists mainly of 3 steps: study the map <fit, 
then solve the PDE (JHD, and show that the generating function of a family (q) solves (1141) if 
and only if (q) solves Smoluchowski's equation (JTJ). The conclusion is then easy to obtain. 

3.2 Inversion of the mapping 

In this section, we study the map 0t, which is useful both for solving theorically the PDE, and 
for obtaining explicit solutions. We will need two preliminary lemmas. 

Lemma 4. Let a > 0, (3, 7 > and K = [0, a] x [0, 0\ x [0, 7] . For (r, s, t) 6 K, denote 

Then for every e > 0, there is a norm ||.|| on M 2 such that 

max \\A(r, s, t)\\ < a + + e, 

(r,s,t)eK 

where we also denote by ||.|| the induced norm on the 2x2 matrices. 

Remark 3. This is a uniform version of the well-known result (see e.g. [25J) which states that 

• For every (square) matrix A and norm ||.||, one has \\A\\ > p(A), where p(A) is the spectral 
radius of A, 

• For every matrix A and e > 0, there is a norm |.| such that \\A\\ < p(A) + e. 
Note indeed that a + \fj3j is the spectral radius of A(a, f3, 7). 

Proof. 1. First assume that (3 and 7 are positive. We can diagonalize A := A(cx, [3, 7). If we 
let a := a, b := y/]3 and c := ypy then 



4 p \" n \\ r 

a — be 



where 



b -b\ „_! 1 / c 6 



c c J ■ ' 26c V-c & 



Now, consider the following norm: for x G M?, let ||x|| = \\P 1 x|| 00 , where ||(xi,X2) 
max(|xi|, \xz\). Then for any 2x2 matrix M, 

M = max - - — — = max— — n = max — — = P MP L. 

An easy computation shows that for (r, s, t) e f^, 



Recall that for a matrix M, 



IMI 



max 



so that, since r > 0, 



P _ M(r,s,*)P| 



6t cs 

r H h — + 

2c 26 



6t cs 
2c ~ 26 



:=F(r,s,t). 



It remains to find the maximum of F on K. First, note that for (r, s, t) G If, 

< F(r,s,t) < F(a,s,t). 

Then, for every (s,t) G [0,(3] x [0,7], we can write t = ps, p > 0. If p < c 2 /6 2 , then 
cs/(26) > 6t/(2c), so that F(a,s,t) = a + cs/6. But s < 6 2 , so F(a,s,t) < a + be = 
a + y/^y. And if p > c 2 /b 2 , then cs/ (26) < 6t/(2c), so that F(ct, s,t) = a + 6t/c. But 
i < c 2 , so F(a, s, t) < a + 6c = a + vT^T- Finally, the maximum of F on F', i.e. the 
maximum of \\A(r, s, t)\\ on F, is a + v7^7- 



2. Assume now that (3 or 7 is zero, say e.g. 7 = 0. Take e > 0, and M > such that 
/?/M < e. Consider the norm = 1 1 Px ||oo, where P is a diagonal matrix with diagonal 
(1, M). For (r, s, 0) G -ft', we have as before 



-1 



\\A(r,s,0)\\ = \\PA(r,s,0)P 



Since s < /3, this shows that ||A(r, s, 0)|| < a + e. 



r s/M 
r 



□ 



We will deal often with real-analytic functions in the remaining of the proofs. For the 
definitions and results on this topic, we refer to [22]. We will show the following result. 



Proposition 1. For t G [0, T c ) and z G [0, 1], define (f> t (., .,z): [0, l] 2 -> IR 2 by 

(j) t (x,y,z) = 

and let K t (z) be the closed subset of [0, l] 2 : K t (z) = <j> t {., ., z) -1 ([0, l] 2 ). Then 



1 + t)x - t-g-to y, (1 + t)y - t-Q^yXi Vi z > 



:i2) 



(i) (f) t (.,.,z) : K t (z) — > [0, l] 2 is a homeomorphism. Denote h t (., ., z) = (h[ 1 \., ., z), h\ A> (., ., z)) 
its inverse. 



,(2), 



(ii) For i — 1,2, (x, y, z, t) 1— ► 6,^ (x, 2) is an analytic function on (0, l) 3 x (0, T c 



Proof, (i) Fix some z G [0, 1] and some t G (0,T C ), and keep the notations of the statement. 
For notational simplicity, we omit the parameter z. Let < t < T c . We first want to 
show that 4> t : K t — > [0, l] 2 is one-to-one and onto. Fix (u, v) G [0, l] 2 and let us check 
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that there is a unique couple (x,y) G [0, l] 2 such that (j) t (x,y) = (u,v). This requirement 
is equivalent to finding a unique fixed point to 



F t (x,y) 



u 



1 + t 1 + t dy 



t dgo, . v 
x,y,z), 



t dg 



1 + t 1 + t dx 



Because of the remark above, F t is a mapping from [0, l] 2 to [0, l] 2 . It remains to check 
that it is contracting. Its differential is 



DF t (x,y) 



1 + t 



/ d 2 g &go\ 

dxdy dy 2 

d 2 go d^gp 

\ dx 2 dxdy ) 



t 



a{x,y,z) P(x,y,z) 
1 +t \lf(x,y,z) a(x,y,z) 



(13) 



Let a = a(l, 1, 1) = (ab, c ), (3 = (3(1, 1, 1) = (b 2 , c ) and 7 = 7(1, 1, 1) = (a 2 , c ). Since 



t < T r , then 

1+t 

there is a norm II. II such that 



a + v^7 + e) < 1 for some small enough e > 0. Hence, by Lemma 0J 



max \\DF t (x, 



< 



1+t' 



a + V/57 + e) < 1, 



so that F t is contracting. Hence it has a unique fixed point. As a consequence, there is a 
unique couple (x, y) G [0, l] 2 such that 4> t (x,y) = (u,v). Moreover, since F t is continuous 
with respect to (u, v) and uniformly contracting in (u, v), then the mapping (u, v) h-> (x, y) 
is continuous, that is h t : [0, l] 2 — * K t is a homeomorphism. 

(ii) For t £ (0, T c ), zq G (0, 1) and (xq, yo) G U to , the matrix Dcj)t (xo, yo, zq) is invertible. Then 
Theorem 2.5.3 in [22] shows that the inverse mapping of 4>t has real-analytic coeeficients, 
i.e. are real-analytic functions on (0, l) 3 x (0,T C ). 

□ 



3.3 Study of the PDE 

The following (non-quasilinear) PDE is a central feature of our discussion 

dgt = dgt dg t 1_ fdg t + y ^9t\ ^ 

dt dx dy 1 + 1 \ dx dy J 

A preliminary result to the proof is the following. Its proof is exactly the same as the one of 
Lemma [TJ 

Lemma 5. Let (ct)te[o,T) be a solution to the system 

^°t(p) = \ ^2p-(p\p') c t(p')ct(p\p) - Y^°t(p) (15) 

for p = (a, b, m) G S , with nonnegative initial conditions. Then for all t G [0, T) and p G S , 
c t (p) > 0. 
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Proposition 2. (i) For every z G [0,1], the PDE (1141) with initial conditions go = go(.,.,z) 
has a unique regular (in the sense of Definition^) solution (t,x,y) i— > gt(x,y,z) defined 
on [0,T C ) x (0,1) 2 . 

(ii) The solution of the PDE is given by 

g t (x, y, z) = g (h t (x, y, z),z)- J^^( h t( x , V, z ), z )^f-( h t( x , V, z ), z ), ( 16 ) 
where h t is defined in Proposition^ 

(iii) We have the alternative expression 

9t{x,y,z) = ^ ^ 2 \x,y,z)+H^(0,y,z^+G t (z) (17) 
in the notations of Theorem^ 

(iv) For every t G [0,T C ), g t has an analytic expansion 

g t (x, y,z) = Y^ b, m)x a y b z m (18) 
for (x, y, z) G [0, l) 3 , where c t (a, b,m) > 0. 

Remark 4. Formula (jT7|) will be useful to compute explicit solutions, since with it, it is enough 
to have the analytic expansion of ht around to obtain the one of gt (whose coefficients are 
precisely the solution to ([1])). Note however that Go m &y be tedious to compute in general, but 
since it is a function of z only, it is relevant only when we wish to compute the concentrations of 
particles with no arms. Nonetheless, their concentrations can be obtained thanks to the sytem 
([T]), since 

d m ~ ^ 

-q(0, 0,m) = Yl OtO-' °> m 'K(0, l,m- m'). (19) 

m=l 

Proof. We will prove the statement in three steps. Fist we will show that a solution has to be 
written as in ( fl6l) . Next that this formula does provide a solution. Proving formula ( flTI) is then 
an easy matter. In all the proof, some z G [0, 1] is fixed. 

1. Let U t = 0t(., ., 2) _1 ((0, l) 2 ), and consider g t a regular solution of (Ti4"j) on [0, T) x (0, l) 2 . 
Fix t £ (0, T) and (x, y) G U to , and let 

{pi{t),P2{t)) := (^{(j)t{x,y,z),z),^(<f> t (x,y,z),z) 

It is easy to see that U t decreases with t, so for t < t , this definition makes sense and we 
can differentiate pi. The regularity assumptions on g t are just those needed to allow the 
use of Schwarz's theorem, and an easy computation shows that, on [0,to], (^1,^2) solves 
a linear differential system with continuous coefficients, whose solution is given by (Q. 
Hence 

^((f) t (x, y, z),z) = ^(x, y, z)/{l+t) ; |p(0t(z, V, z ), z ) = ^ > z )l ( 1+t ) ( 20 ) 

for all (x, y) G U t . Then, it is easy to check that for all (x, y) G U t , ffTUj) and (TTTT) hold. 
Replacing (x, y) by h t (x, y, z) (recall h t : (0, l) 2 — > U t is the right-inverse of (fi t ), we finally 
obtain fTi75]) . This shows that the PDE has at most one solution. 
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2. The existence of a solution is now straightforward. Let gt be defined as in (1161) . Because 
of the regularity of h t and of go, g has the required regularity properties. It then suffices 
to show that it is actually a solution. To this end, let us first compute 

(pi(t)Mt)) ■= (t^)'§^)) 

for some fixed t G [0,T C ) and (x,y) G U t . By differentiating gt(4>t) with respect to x and 
y, it is easy to see that it solves a linear system, which has, before T c , a unique solution, 
given by equation ([20]) . To conclude, we may differentiate gt(4>t) in two different ways: one 
using ( fl6l) and ( !20l) . The other with the chain rule. Compounding by h t in the obtained 
equality readily shows that g t solves the PDE ( fl~4"|) for t G [0,T C ), (x, y) G (0, l) 2 . 

3. The formula (T5]) is easy to obtain, by differentiating g (h t (x, y, z), z) with respect to x, y 
and z, and using the fact that 

^9o, u \ r(2) ®9o f i \ t(i) 

— {ht, z) = h\ > ; -j-(ht, z) = h\ ' . 
ox ay 

in the notations of Theorem [TJ 

4. To prove the last point, consider t G [0, T c ). (f> t is well-defined and analytic (in (t, x, y, z)) 
in a neighbourhood of (t , 0,0,0), and _D0 to (O,O,O) is invertible. So, by theorem 2.5. in 
[22] . /it is analytic near (t , 0, 0, 0), hence so is <ft = go(h t ). So we may write 

g t (x)= c t (a,b,m)x a y b z m (21) 

for (t,x,y,z) in a neighbourhood of (to, 0,0,0) and infinely differentable (even analytic) 
c t . By analytic continuation, the (c t ) are uniquely defined, so we can let 

E = {t £ [0,T c ),Vp G 5* ct(p)>0}. 

By continuity, E is a closed set containing 0. On the other hand, (121]) holds for (t, x, y, z) 
in a neighbourhood of (t , 0, 0, 0), so for t G E, there is a e > such that ff2T]) holds for 
t G (to — e, + e) and (x, y, z) G (— e, e) 3 . In particular, since g t solves the PDE ( f!4l) . it 
is easy to see, using a Cauchy product and identifying the coefficients, that c t solves ( |T5i) 
for t G (to — e? to + e) • So by Lemma (t — e,t + e) C E. So is open and E = [0, T c ) . 
Finally, recall from Proposition [1] that h t , and so gt, are analytic on [0, l) 3 . But we have 
just shown that gt has an analytic expansion around with nonnegative coefficients. So 
(see e.g. the proof of Berstein's theorem in [22]), this expression actually holds on [0, l) 3 . 

□ 

3.4 Equivalence between the system and the PDE 

Smoluchowski's equation is solved thanks to the PDE (1141) . 

Proposition 3. (i) Let (ct)t£[o,T) be a solution to Smoluchowski's equation, and let gt(x,y, z) := 
(c t ,x a y b z m ) be its generating function. Then for all z G [0,1], (t,x,y) i— > g t (x,y,z) is a 
regular solution to the PDE (fT4"|) on [0,T AT^) x (0, 1) 2 ; with initial conditions g$(., .,z). 
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(ii) Conversely, let {c t (p)) P es, t G [0, T) be a family of differentiable functions. Let g t (x,y, z) 
be its generating function and assume it is defined for t G [0,T), (x,y) G (0, l) 2 and 
z G [0, 1] . Assume that for every z G [0, 1], g t (., .,z) is a regular solution to the PDE ffT4l 
with initial conditions g (., ., z). Then 

• For allpe S and t G [0, T), c t {p) > 0, 

• (ct) is a solution to Smoluchowski's equation fort G [T AT C ), with initial conditions 

Co- 

Remark 5. An important feature of this result is that the PDE (fT4"l) and the system ([1]) are 
equivalent only before the critical time (T c or T^). This fact is crucial when we study the 
microscopic model. We indeed obtain a family of coefficients whose generating function solves 
the PDE (on [0, +oo)), but we cannot ensure that they solve Smoluchowski's equation after T c 
(actually, we believe that they do not). 

Proof of Proposition^ (i) First note that g is regular according to Remark [2J If one takes 
f(a, b, m) = x a y b z m in for some fixed (x, y, z) G (0, l) 2 x [0, 1], then one gets 

d ( v dg t dg t dg t dg t 

-T79t{x, y, z) = — A t y- B t x — . 

at ox ay ay ox 

Recall from Lemma[3]that when t < T^, A t = B t = 1/(1 + 1). Replacing in the equation 
above shows that g t solves (fT4"l) for (x, y) G (0, l) 2 and < t < T A Y^. 

(ii) As in the fourth part of the proof of Proposition (2], we see that the (c t (p)) solve ffl"5l) . 
and hence that they are nonnegative. By uniqueness of a solution to the PDE f|14l) . 
for t G [0, T A T c ), these coefficients are those obtained in f[T8l) . Now, let t < T c A T, 
U t = 0i _1 (., ., 1)((0, l) 2 ), K t = K t (l) = 4>~\., ., 1)([0, l] 2 ), and recall from that since 
g is a regular solution to (JHJ), then for all (x, y) G U t 

^(»,v,i),D = ^(»,»,i) 1 



ax ax 1 + 1 

what we can write 

£ a Ci (a,6,m)#(x,y,l)>| 2) (x,y,l) b = £ ac (a, 6, m)xV x (22) 

(a,b,m)gS (a,b,m)£S 

Note now that since £ < T c , then </> t (., ., 1) : K t — > [0, l] 2 is a homeomorphism, so C/ t = K t . 
Since (1,1) G K t , we can pass to the limit in the equality above when (x, y) — > (1,1). 
Using monotone convergence and the continuity of fa, we obtain 

£ ac t (a,b,m)= £ ac (a, b, m) x = 

(a,b,m)£S (a,b,m)£S 

The same reasoning shows that (6, q) = 1/(1 +£) for £ < T c . Hence, we may re- write (1151) 
before T c by substituting 

_^_ = a £ b'ct(a',b',m') ; ^— = 6 £ a ct(a , b' , m') , 

(a',b',m')eS (a',b>,m')eS 

which shows that (c t ) solves Smoluchowski's equation before T c . 

□ 
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3.5 Existence and uniqueness of a solution 

With these results, proving Theorem [I] is now an easy matter. 

Proof of theoremUi 1. Let us first prove that (a 2 + b 2 , c t ) is finite before T c and tends to +oo 
when t — > T c . So take (c t (a, b, Tn))te[o,T) a solution to the system (TjQ), and g t its generating 
function. Since (a 2 + b 2 , c t ) < +oo in a neighbourhood of 0, then we have 



&9t 
dx 2 



(1,1,1) = (c t ,a -a) 



as long as (a 2 + b 2 , c t ) < +oo. Note that by Lemma [3] (c t , a) is bounded by 1, so (c t , a 

dx 

[ 1201) with respect to x and y and having (x, y) tend to (1, 1), we obtain that 



explodes if and only if -?^-(l, 1,1) explodes. Let us compute the latter. Differentiating 



1 + t — ta —tj \ f a \ 1 ( 7 



l + £ — toy \c J 1 + t V a 



where 



« = ^ = <aW ; P = W = (b -6, co> ; 7 = ^ = (a - «, <*>. 

and 

■ =§£(1,1,1) and »=^(1,1.D. 

Hence 

a=d S {W)= S a(a-l)c t (a,&,m)= ^ 

(a,fe,m)G5 V 7 

This expression is valid as long as t < T c , since the determinant of the matrix is then 
nonzero. In the same way, we also have 

^(1,1,1)= £ Kb-lMaAm)^ f 

If 7 or (5 is nonzero, then (c t , a 2 + b 2 ) — > +oo when £ — > T c . If 7 = /? = 0, then (a 2 + 5 2 , c t ) 
remains finite, but this condition also imposes that M — 1, and so T c = +00. 

2. Uniqueness is now easy to obtain: assume (c^) and (c$ ) solve the system CD on [0,T), 
T < T c , with initial conditions cq. Let gl and # t be their generating functions. Since 
Too = T c and T < T c , then by Proposition [3j for every z G [0, 1], they are regular solutions 
to the PDE (fl4j) on [0, T) x (0, l) 2 , with initial conditions g (., .,z). But by Proposition 
[2] there is a unique regular solution to the PDE on [0, T c ), so g^' = on [0, T), so that 

( C | 1} ) = ( C | 2) ). 

3. The existence is given by Proposition [2], (iv), and Proposition El (ii). 
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4. Let us finally prove that the total mass is conserved. Consider ^(x, y, z) = {4>t{x, y, z), z), 
U' t = ■0 t " 1 ((O, l) 3 ) and K[ = ^ t _1 ([0, l] 3 ). For (x, y, z) e U' t , we can differentiate g t (ipt{x, y, z)) 
with respect to z, and using fl20|) . we obtain 



^(ijj t (x,y,z)) = ^(il>t(x,y,z)). 

Now, ipt is a homeomorphism from K[ to [0, l] 3 , so U' t = K t . But (1, 1, 1) e so we 
may pass to the limit when (x, y, z) — > (1, 1, 1) in the equality above, to obtain 

^(1 1 1) = ?*(1 1 1) 
what precisely means (m,c t ) = (m, Co). 

□ 



4 Explicit formulae 

We give in this short section some explicit solutions, without giving the full details of the 
computations. We will always assume that at time 0, there are only particles of size 1 in the 
medium. So given a (finite) measure \i on N x N, we assume 

co(a,b,m) = ii(a,b)l {m=1} 

and as usual (a, Co) = (b,co) = 1. To obtain the solutions, we need to invert <p t , what can be 
done using the (two- variable) Lagrange inversion formula (a statement is given by Good [12]). 
But it is much more involved than the one-dimensional formula, and the expressions it would 
provide can hardly be called explicit. Let us however study three easy cases. Only the last one 
requires the two-variable formula. 



4.1 Particles with one female arm 

The first case is when each particle has exactly one female arm, and a number of male arms 
distributed according to a measure H\. So 



H(a,b) 



if b ^ 1 

/ii(a) if 6=1 



and we will assume that Aq = Bq = 1, i.e. \i\ is a probability measure with unit mean. In this 
case, we obtain, for every a, b > and m > 1 

if b / 1 

ct{a,b,m) = { t m - x 1 fm + a — 1\ 

* ' Vi i m + «.-!) if b = 1 



(1 + t) m+a m\ a 

In particular there exists only particles with one female arm, what is physically obvious. More- 
over, the concentration c t (a, 1, m) is exactly the concentration of particles with a arms and mass 
m obtained in the "oriented model" of [2], with initial distribution fii. This is also natural, 
since in this case, (a, 1, m)- and (a', 1, m')-particles indeed coagulate with rate a + a', which is 
the rate of the oriented model. Note also that T c = +oo, like in the oriented model. 
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4.2 Arms with uniform random genders 

In this model, the total number of arms of a particle is chosen according to a measure /ii, then 
each arm is given a gender independently, with probability 1/2. That is, we let 



fj,(a, b) = fii(a + b) 



a + b\ 1 



b 2 



We will assume that /ix has mean 2, so that A = B = 1. Let = (J + l)/i(j + 1). Then we 
obtain, for (a, b) ^ (0, 0) 

1 (rn^a + b-2)! my™ 

Q(a,o,m) = -- . , , , — — — (m + a + b — 2) 

V ' ' 7 2(1 +t)a+6+m-i m \ a \ h \ \ 2 J v ; 



and 

c t (0,0,m) = — — -(— ) (m-2) 

V ' ' ; m(m-l)2(l + l/t) m - 1 \2J V ; 

provided ^i(O) > 0. If T c = +oo, this condition means that v\ ^ 8\. In particular, one easily 
checks that 

^2 c t (a,b,m) = 2cf m (fc,m) 

a+b=k 

where c s t ym (k,m) is the concentration of particles with k arms and mass m in the symmetric 
model of [2J, with initial arm distribution /ii/2. The factor 2 comes from the normalisation: in 
our model, the total concentration of arms in the medium is 2, when it is 1 in the symmetric 
model. It is also worth stressing the stronger fact that for a, b > 0, we have 

ii \ 1 fa + b\ 1 
ct(a, b, m) = - I b J ^ct{a + b,m). 

Hence, at any given time, the distribution of the number of male (or female) arms is still 
binomial. So, if at some time we chose to reassign to each arm a gender uniformly and inde- 
pendently, and let the system evolve on from this state, no difference would be observed. Or 
we could watch a system evolve like the symmetric model starting from an arm distribution 
/ii/2, and then at some time give the arms a gender uniformly at random and independently. 
The evolution afterwards will be the evolution of the sexed model with initial arm distribution 
fi. Note as before that the critical time is the same than in the symmetric model with initial 
distribution yUi/2. 

More generally, consider initial concentrations such that for all a, b > /i(a, b) = fi(b, a), and 
the solution {c t )t^[o,T c ) to Smoluchowski's equation ([T|). Then it is easy to check (by uniqueness) 
that for all t e [0, T c ) c t (a, b, m) = c t (b, a, m), and that, if we denote 



k t {l,m) := 22 cti^^m), 



a+b=l 



then (k t ) is governed (up to a factor 1/2) by the symmetric Smoluchowski equation of [2]. 
Hence, in this case too, k t (l,m) = 2c s t ym (l,m). 
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4.3 Particles with one gender 



Let us finally consider the more intricate case where at time 0, the arms of each particle have all 
the same gender. This is motivated by the idea of ionic bonds: a particle with only male (resp. 
female) arms can be considered as a cation (resp. an anion), and cations can only bond with 
anions. Hence, consider, for i — 1,2, fii two measures with mean 1 such that ^\(0) = /i2(0), 
and take 

{/ii(a) if 6 = 

/i 2 (6) if a = 

else 

and Vi(j) = (J +1) fj,i(j +1) . The two-variable Lagrange inversion formula gives, for (a, b) ^ (0, 0) 

, , v t™- 1 ^(m-k + b-l)\(k + a-l)\ *{ m -k), u , ^ * k( , 1|M 

c t (a,b,m) = — -r — ——7 — r / t , , ,, , , T , iV '[k+a—l)v 2 (m—k — l + b). 

K ' (1 + +a+ ^ (m - k)\k\a\b\ 1 v 7 2 v ' 

If we also let 

m— 1 



z/i ov 2 {m) = (m - 1) ^ T^i (m - r (/ 2 



m — k 

k=l 



then, for m > 2, 



Ct(0,0,m) 



1 1 



m - 1 (1 + 



— rr -i/i o Uo(m), 

l/t) m 



provided i^i(0)i/ 2 (0) > (which means, if T c = +oo, that z/ x and v 2 are not 5i). In particular, 
we see that if T c = +oo and v\, v 2 ^ 6%, then for m > 2, 



ct(a,b,m) 




if (o,6)^(0,0) 
= 0. 



Hence all the arms are used to coagulate. Chemically, this means that there are no more ions 
in the medium. The limiting distribution of the sizes is given by (m — l) _1 ^i o ^(m). We 
will generalize this fact in the following section, and give a probabilistic interpretation of the 
measure Ui o vi- Also, if is the mean of U{, then M = \fMiM 2 and T = 1/(M — 1), or +oo 
if M < 1. If \i\ = \l<i = //, then the critical time is the same as in the symmetric model with 
initial distribution //. 



5 Limiting concentrations and Galton- Watson processes 
5.1 Convergence of the concentrations 

In this section, we will study the limiting concentrations. Similarly to what happens in the 
oriented and symmetric model of [2], we expect the concentrations to converge when the time 
tends to +oo, whenever gelation does not occur. Physically, this would mean that the system 
converges to a terminal state where all arms have been used (otherwise, further coagulations 
"should" occur). This is actually true, and this is an easy consequence of the preceeding results. 

Corollary 1. Assume T c = +oo, and let (ct)t>o be the solution to Smoluchowski's equation (0Q). 
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(i) When t — » +00, there exists limiting concentrations c^m) such that 

ct(a,b,m) -> 000(771)1^=6=0} 

m ^(N x N x N*). 

(ii) For 2 G [0,1), the generating function g^z) of (c x (m)) m >\ is the antiderivative vanishing 

at of 



where (hx,hoo) is characterised by 

h£\z) = f^(h£\z),hg\z),z)) 



#/ V 7 (23) 



Proof, (i) Since = T c = +00, then ([2]) holds for all t > 0, so 



Ec t (a,b,m) — > 0. 
/ — s--t-rm 

■6)^(0, 

Then, using (TT91) . we get for all t > 



f — + oo 

(0,6)^(0,0) 



m-1 „ t 

Cj (0, 0, m) = cq(0, 0, m) + / ] / c a (l, 0, m')c s (0, l,m — m') ds. 

m'=l ^ 

But the integrand is bounded by A S B S = 1/(1 + s) 2 . Hence the integral has a finite limit 
when t — > +00, and so does q(0, 0,m). Finally, (m,c t ) is bounded by Theorem [lj and 
(m, Coo) < +00 by Fatou's lemma, so the Cauchy-Schwarz inequality shows that 



E|c t (0,0,m)-Coo(m)| — > 
t— >+oo 



m>l 

and the result follows, 
(ii) By ^-convergence, we have 

000(2) = lim g t (0,0,z), 

t— >+oo 

so, using ([2} and the fact that i/^ and ^ 2 ' ) are bounded by 1, 

Ooo(^) = lim G t (z). 

It just remains to check that h[ (0, 0, z) and /i^(0, 0, z) do have a limit when t — > +00. 
From their definition ([3]), they have the same limit (if any) as k\ (z) := h t (0, 0, 2;) and 
fcj ; (^) := ; (0, 0, z). Now fix Z e [0, 1), and consider k\ and fcj as (continuous) maps 
on [0, Z\. But h t is the right-inverse of t , so 

^(tf'W, *f W, *) = if W ; f^M, *P>M, ») = tf'M (24) 
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and, since (a 2 + 6 2 ,co) < +00, —— (x, y, z) has a bounded differential on [0, ll 2 x \0,Z]. 

ox 

Hence kf\ and for the same reason k[ , are Lipschitz-continuous on [0, Z], with a constant 
independent of t. Ascoli's theorem thus shows that the families {k\ ) and {k\ ), £ > 
lie in a compact set (for the uniform topology on [0, Z}). So the family {k[ , k^) lies in 
a compact set, and passing to the limit in fpM]) shows that any of its limit points solves 
( 1231) . But since T c = +00, the application 

f \ ( dg °( \ dg ° ( \\ 
[x,y) i-> I — {x,y,z),— \x,y,z) I 

is contracting for every z G [0, Z]. So there is a unique solution to (|23l) . and (/c^,/^) 
converges to this solution. 

□ 



5.2 Connection with two- type Galton- Watson processes 

In [2], Bertoin shows that for monodisperse initial conditions (i.e. co(a,m) = /x(a)l{ m= o} 
for some measure /i) and when gelation does not occur, the limiting concentrations can be 
described in terms of Galton- Watson processes. The same kind of analogy is observed in our 
case. Precisely, consider a Galton- Watson tree with two genders, constructed as follows. We 
start from a male or a female ancestor. It gives birth to a number a of male children, and 
a number b of female children, where (a, b) is distributed according to a law /j, m (a,b) if the 
ancestor is a male, /i/(a, b) if it is a female. Then each child gives birth to a certain number 
of children, distributed according to \i m or /j,f, depending on his gender, and so on. Consider 
Tf(fjL m ,jj,f) (resp. T m (p m , }if)) the total population of such a Galton- Watson process, starting 
from a female (resp male) ancestor. Let for r G [0, 1) 

g f (r) = E(r T f) ; g m (r) = E(r T ™) 

their generating functions. It is an easy exercise to check that they solve the following system, 
where 0/ (resp. <f> m ) is the generating function of [if (resp. /i m ) 

f 9m(r) = r(f> m (g m (r),g f (r)) 
I 9f(r) = r(f>f(g m (r),g f (r)). 

Besides, if T(/i m ,/i/) is the size of a Galton- Watson tree started from a male and a female 
ancestor (each tree growing independently), then 

P(T(> m) n f ) = m) = [z m ]g m (r)g f (r), (26) 

where [z m ]h(z) is the coefficient of z m in the expansion around of an analytic function h. 

Now, let us go back to our study. Assume monodisperse initial conditions, i.e. there is a 
finite measure /x on N x N such that co(a, b, m) = l{ m=0 }/i(a, b), and assume T c = +00. We will 
use the same notations as in the previous section. Using (1231) . we obtain 

9M = -h£(z)h<£(z), 
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so that for n > 2 



[z n \9oo(z) = ^[z n ]hU(z)h£\z). (27) 

v m (a, b) = (b + l)/x(a, 6 + 1) the probability measure with generating fuction (f) r 
ff(a, b) = (a + l)/x(a + 1, b) the probability measure with generating fuction <pf 



Let 

. = dgo 

dy ' 

= <%o 

dx 

(they are probability measures because of the assumption (a, c ) = (b, c ) = 1), and consider the 
two-type Galton- Watson process as above with these reproduction laws. Because of (|25|) and 
(1231) . (g m ,9f) an d , solve the same equation, which has a unique solution by Corollary 
HI so h ( oo ] = g m , h£ ] = g f . Hence by ([26]) and (J27J) 

P(T(i/ TO) i/,) = m) = [^/^(z)/^) = (™ - l)[* m ]<7oo(*) = (m - l) Coo (m). 

Finally, let us call a measure /x degenerate if /x = 5(1,1) or /x = -(5( 2 ,o) + ^(0,2)), or /x(a, 6) = 

for a 7^ 1, or /x(a, 6) = for 6 ^ 1. We let the reader check (using e.g. Theorem 10.1 in [T7] ) 
that under the assumptions T c = +00, and ruling out the degenerate cases, T m (v m ,Vf) and 
Tf(fjL m ,fif) are finite a.s., and that the process is supercritical if T c < +00. 

Corollary 2. The limiting concentrations verify for m > 2 

P(T(i/ m , i//) = m) = (m - l)coo(m). 

Moreover, if /x zs noi degenerate, then T{v mi Vf) < +00 a.s. 

Hence, the law z/i o z/ 2 defined in Section H~3l is the law of the total population of a two-type 
Galton Watson process started from one male and one female ancestors, where the males give 
birth to females according to the law U2, and the females give birth to males according to the 
law v%. In particular, if v\ = = v, then vov is the distribution of the size of a Galton- Watson 
tree with reproduction law v and starting from two ancestors. So we get (what is not obvious 
from the formula for o), that for m > 2 

2 

v o u(m) = —v* m (m - 2). 
m 

This corollary answers another question about gelation. By Theorem HI the total mass (m, q) 
is conserved as time passes, so gelation does not occur before T c . But if T c = +00, it may occur 
at infinity: some mass may be lost then. For monodisperse initial conditions, Corollary [2] proves 
that this cannot happen, except in the degenerate cases. Denote indeed C t = (c t , 1). Because of 
the ^-convergence in Corollary [T], Ct — > Yl c oo(^), and Equation (jSJ) yields C t = Cq — t/(l+i). 
Hence 

C - 1 = ) j Coo(m). 

m>l 

So, whenever /x is not degenerated, Corollary [2] gives 

} u mcoo(m) = ^ (("i - l)coo("i) + Coo(m)) = P(T(/i m ,/i/) < +00) + ^Coo(m) 

{a,b,m)£S (a,b,m)eS 

=1 + Yl c °°( m ) = 1 + Co - 1 = C , 

which is precisely the total mass at time 0. For the degenerate cases, we can get explicit 
expressions for the concentrations, and these show that the mass at infinity is 0. 
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6 Microscopic model 

6.1 Notations and preliminary results 

The goal of this section is to construct a sequence of random processes modeling the coagulation 
of particles with male and female arms. We will start with n particles (and then let n — > +00). 
Let us first set some notations. 

• Recall S = N x N x N*. 

• [0,n] = {(),..., n}. 

• M > is a fixed real number. The number of arms and the total mass are assumed to 
grow at most like Mn (see the definition of E n ). 

• For p = (a, b, m) G S and p' = (a', b', w!) G S, we will denote p.p' = a'b + ab' the rate of 
coagulation, and pop' = (a + a' — 1,6 + 6' — l,m+ m') the type of the particle resulting 
from such a coagulation. 

• The sequence of the number of p-particles is an element of 

E n = I N G [0,n] s , (a + b + m)N(a,b,m) < Mn 

^ (a,b,m)eS 

which is a finite set. 

• —E n is a. subset of 
n 

E = <j C G [0, l] 5 , (a + b + m)C(a,b,m) < M 

(a,b,m)£S 

An element of E represents the sequence of concentrations of p-particles. E is a metric 
space endowed with the distance 

d(C^ l \ c^) = \c {1) (p) - c (2) (p)| • 

pes 



• We will call C-convergence the compact convergence (i.e. uniform convergence on every 
compact set) for functions from M + to E. 

• D([0, +00), H) is the space of cadlag functions from [0, +00) to a metric space (H,d), 
endowed with the Skorokhod distance. We will call S-convergence the convergence for 
Skorokhod's distance. For the basic facts about Skorokhod distance for functions with 
value in a (complete separable) metric space, see [7]. 

It is easy to check the following result. 

Lemma 6. (E, d) is a compact metric space. In particular, it is a Polish space. 
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6.2 Model 



Let us now introduce the model. Informally, we consider a finite number n of particles with 
integer mass, and assume that at time 0, the total mass of the system plus the total number of 
arms is less than Mn. Then, each pair formed of a p-particle and of a p'-particle may coagulate 
with rate \p.p' ', independently of the other pairs, to form a p o p'-particle, that is, the time one 
has to wait to see them coagulate is exponential with parameter \p.p' ■ In other words, assume 
the system in in the state rj at a given time, that is i] e E n and rj(p) is the number of p-particles. 
There are r](p)r](p') (or r)(p)(r)(p) — 1) if p — p') pairs formed of a p-particle and of a p'-particle. 
Let 

\p-P't]{p)t]{p') if p 7^ p' 

\p-pr){p){r)(p) — 1) if p — p' . 



\{P,P') 



We set independently on each couple (p,p') an exponential clock with parameter X v {p,p') (an 
exponential random variable with parameter is assumed to be a.s. infinite). There is a.s. one 
and only clock which rings first. If it is the clock on the couple (p,p'), then the system jumps 
to the state 77 + A p>p i where 

A p,p'(p) = A p,p'(p') = -1 if P^P' 
A p,p'(p) = - 2 if P = P' 

\, P '(P°P') = +1- 

Then restart the construction afresh from the new state. Note that only finitely many rj(p) are 
nonzero, so the first jump occurs after an exponential time with parameter 

\ v = \(p,p') < +00. 
p,p'es 

We will consider the Markov chain constructed according to this rule. That is, we fix for every 
n > 1 

• An element Xq^ of E n , which is the initial number of particles. 

• A pure-jump Markov process IW on E n , defined on some probability space (f2 n , A n , P n ), 
starting from X^ n \ and with generator 

Gf(v)= E (f(v + \, P ')-f(v))\(p,p') 

(py)e5 2 

for every bounded function / : E n — > R. The construction of such a process is obvious 
since E„ is finite. 



The rescaled and time-changed process 



n 



Note that is a pure-jump Markov process on ^E n C E, starting from = Xq^/ti, and 
with generator 



(p,p')GS 2 
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where 

Ai n W) = ^A n >,p'). 

The law P n of the process is a probability measure on ©([0, +00), E). We will prove that 
the sequence (P n ) is tight, and that for every limit point P, and almost every process (C f ) 
with law P, (C t ) solves some system, which is Smoluchowski's equation (pQ) before the critical 
time. Because of the uniqueness of such a solution, this will show that (P n ) itself converges 
to the solution of Smoluchowski's equation before the critical time. The proof of tightness is 
analoguous to the one in [18], up to some slight modifications. 



6.3 Tightness 

Lemma 7. The sequence (P n ) n >o is tight. 

Proof. We will use the classical tightness criterion stated in [19] . page 34, or in [7], Theorem 
7.2. For t > 0, let P t {n) be the law of C^ n \ which is a probability measure on E. Since E is 
compact by Lemma El the tightness of the sequence (P t ) n >o is obvious. 

Now, C^ n ' is a pure-jump procees on -E n C E, with generator G^ n >. Hence, when the 
process is in the state 77, then the time before the next jump is exponential with parameter 

A J n) := Y x ?(pip') = 2 ( Y n P-p'v(p)v(p') - ^p-vnip) 

P,p'es \p,p'es pgS 

and, since rj e E, A$ n) < M 2 n := cn. Now take iV > 0, (3 > 0, e > 0, and let 5 > such that 
N = 51 for some I £ N*, and 3c5e//3 < 1. Define now 

w N (Y,5):= inf max sup d(Y t ,Y s ) 

wElls U£n ti<s<t<U +1 

lis being the set of all subdivisions = t < ti < ■ ■ ■ < t n = N of [0, N] such that t i+ i — ti > 5 
for all i. Consider the partition t = < ti = 5 <■■■< ti = \ = N of [0, N}. Let 
Z t := sup rf(d n) , Cl n) ) for < % < I - 1. Then 

ti<S<t<t i + 1 

F n (w N (C (n) ,5) >0)<F n ( max Z i > ) < I max P n (Z; > 

yO<i<i — 1 / 0<i<(— 1 

But the size of a jump, that is d(C^_ , C^), is 3/n. Hence, if Zi > (3, then the process has 
jumped more than k := \/3n/3] := \c'n\ times between ti and ti+i (where \x] is the first integer 
strictly greater than x). If Sk if the time of the k-ih jump, the Markov property tell us that 

But Sk is the sum of h independent exponential random variables, with parameter smaller than 
cn. So, if S' k is the sum of k independent exponential random variables with parameter an (on 
a probability space (Q, A, P)), then is stochastically dominated by S' k , that is 

Pn(5 fc <5)<P(^<5). 

To conclude, note that the last term is the probability that a Poisson process with parameter 
cn jumps more than k times on [0,6], and Stirling's formula shows that this tends to zero for 
(c5e/c') < 1. □ 
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6.4 Convergence 

In this section, we prove the convergence of to a process solving the system ([15]) . and 
deduce that it solves Smoluchowski's equation ([I]) before T c . 

Proposition 4. Assume that the following convergences in distribution hold 

• For every p G S , Cq^Qo) — > Co(p) for some non random Cq(p) > 0, 

• aCp (a, 6, m) — > 1, 

(a,b,m)eS 

. £ bC^ n \a,b,m)^l. 

(a,b,m)eS 

Let P be a limit point of (P n ), and let (c t ) be a process with law P. Then a.s. (c t ) solves the 
system ( flBI) . with initial conditions (co). 

In the following proofs, we take a subsequence of (C^) which converges in law to some 
possibly random c G HJ([0, +oo), E). For notational simplicity, we will assume that (C^) itself 
tends to c. Since E is compact, it is separable, and hence, so is D([0, +oo), E). Skorokhod's 
representation theorem (cf e.g. j7]) now allows us to assume that the C^ n > are defined on the 
same probability space (ft^P), that C {n) -> c a.s. (that is, for almost every lu G fi, the 
function C^ n \cu) tends to c(tu) for Skorokhod's distance), and that in the statement, there is 
a.s. convergence. We will also constantly use the fact that for every bounded Borel function 
/ : E — > R, the processes 

c (n) _ c (n) _ f* G (n) /(C (n) ) ds . = j^n) (2g) 
JO 

and 

(il4 n) ) 2 - (G(">(/ 2 )(CW) - 2/(Cf ))G^/(Cf ))) ds (29) 
are martingales. Note also that if / : E — > R is "linear" , then for all G 

p,p'g5 pes 

and that 

(G ,(n) (/2) _ 2/G H /) fo) = i- ( £ / (A^^^pVp' - i J2f ( A p,p) 2 V(P)P-P ) (31) 

n \p,p'eS U pes J 

We will also need the following convergence result. 

Lemma 8. Let 

= aCi n) (a,b,m) and = ^ &C s (n) (a, 6, m). 

(a,b,m)eS (a,b,m)eS 

Then (A (n) ) and (B^) C-converge a.s. to t \-> 1/(1 + t). 
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Proof. Obviously, we cannot pass to the limit immediately in these expressions. So consider 
the maps from E to R: C h (a, C) and C \— > (b,C), which are measurable and bounded (by 
M). By ([28]), there are martingales M A <^ and M B '( n ) such that 

At ] = 4 n) - f ds + M/' (n) ; B\ n) = B^ - f A^B™ ds + M? ' {n) . (32) 

Jo Jo 

Now, ( l29l) and (I3TT) show that the quadratic variation of M A '^ verifies 

MH 



> M Mn)\<l f A B ds < 



n 



By Doob's inequality, 



Ell sup Mf' (n) 1)^0 



^0<t<T 

for all T > 0. Hence there is a subsequence of (M A '( n )) which C-converges a.s. to 0. In par- 
ticular, it S-converges. For notational simplicity, we will assume that (M A ''™') itself converges. 
For the same reason, we may assume that (M B ^ n ^) also S-converges a.s. to 0. 

Now, note that the proof of Lemma [7] still works for (A^), since the size of its jumps is 
bounded by 1/n. So (iW) is tight, and by Prokhorov's theorem, this means that for almost 
every cu E f2, (A^^lu)),,-^ lies in a compact of D([0, +00), E) (actually, this is a consequence 
of the proof of Lemma 9, and we do not need this implication of Prokhorov's theorem). The 
same works for (B^), so we can find fi' C with P(fi') = 1, such that for every oj E Cl', 
C n (u) -> c(w) (for Skorokhod's distance), M A ^ n \uo) -> and M B '("V) -> (compactly), and 
(vl^ n ^(u;)) n >o and {B^ n \u)) n >Q lie in a compact. 

Next, fix u E Q', and let us find the limit of (A^(u),B^ n '(u>)) (for the product topology 
- which is not Skorokhod's topology on ©([0, +00), E 2 )). Since it lies in a compact set, it is 
enough to show that it has only one limit point. So assume (A^ n \u), B^ n \uj)) converges to 
some (A, B). Then {A\ n \u)) converges to A t for every t E K } the set of continuity points of 
A. But A is cadlag, so it has only countably many points of discontinuity, and hence K c has 
Lebesgue-measure 0. Hence (A^(u)) converges to A Lebesgue-a.s., and ditto for (BW(w)). 
Also, and (B^(u)) are bounded by M, so using dominated convergence in (1321) and 

recalling that Aq(oj) and Bq{uS) — > 1 by assumption, we obtain 

A t = 1 - / A S B S ds ■ Bt = 1 - / A S B S ds. 
Jo Jo 



Hence 



A t = B t 



l + t 



Finally there is only one limit point to (A^ (u) , B^ (u)) . So and (B<")(w)) both 

S-converge to 1 1— ► 1/(1 + 1), and, since this function is continuous, they C-converge. □ 

Remark 6. As pointed out in the proof, the convergence of A™ and B^ to the actual number 
of arms 

A t = ac t (a,b,m) and B t = bc t (a,b,m) 

(a,b,m)eS (a,b,m)eS 



26 



is not obvious. There is no such problem for a strictly sublinear coagulation rate (as in Jeon's 
proof [IB])- In our (linear) case, we prove below that this convergence holds before the critical 
time (we also refer to Norris [21] for general sublinear rates in a model with no arms). In fact, 
if there is a solution (c t ) to ([I]) defined after T c , we believe that A^ 1 and do not converge 
to A t and B t after T c (and that this number of arms is then stricly lesser than 1/(1 + 1)). This 
would suggest that T c is actually a gelation time: some of the arms are lost in a "gel" (a particle 
with an infinite mass and infinitely many arms). 

Proof of Proposition^ 1. Take some p = (a ,b ,m ) G S, and let for C G D([0, +oo), E), 
f(C) = C(p ). According to (128]), 



C[ n \p,) - C ( n \p ) - / GMf(Ct*) ds := Mf' (n) (33) 



is a martingale. Note also that for p,p' G S, /(A p y) is 0, except if p or p' or p o is p - 
Hence, it is easy to check using ( 1281) that 



G (n) /(c (n) ) = _ J- Ci n \p)Ci n \p ) P0 .p P-(P0\p)Ci n \ P )Ci n \p \p) 

(34) 



n 

pes 

The last term is due to the difference between X v {p,p') when p ^ p' and when p = p'. In 
any case, it tends to uniformly on R + and uniformly in p . 

2. Let us now study the martingale term. By Doob's inequality, we have for every T > 

2 s 



E f ^ sup Mf°' (n) ^ j < 4E ((M£°' (n) )^ 



and by ( |29]) . this last term is 



T 



E / (G (n) / 2 -2/G'( n) /)(Cf ) )ds 



But by (HSU), and since / (A PtP ,) < M for all p,p' G S, then (G( n )/ 2 - 2/G( n )/)(d n) ) < 
M 4 /n, so that 

E^sup Mf' (n) ^ -.0. 

Hence, there is a subsequence of (M po,(n )) which a.s. converges to uniformly on M + . For 
notational simplicity, we will now assume that (jkP '( n )) itself C -converges to 0. Using 
the diagonal method, we may as well assume that (M Po,(ra )) C-converges to for every 
Po G S. 

3. We have already seen in the proof of Lemma[7]that d{C[ n \ C._ ) < 3/n a.s. By continuity 
of I h sup sg [ t ] d(X s -,X s ) (cf [7J), this ensures that c is almost surely continuous , so 
actually C -converges to c. From the definition of d, it is also obvious that C <yn \p) 
C-converges to c(p) for every p G S. 
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4. With these results, we may now pass to the limit in (1331 and (1341) . Write (1341) in the form 

G (n) /(c ,(n) ) = _ an{s) + + 

Equation ( 1331) shows that 

Cl n) (po) = C {n \po)- [ a n (s)ds + l [ f3 n (s)ds + f e n {s) ds + M t (n) . (35) 

Jo ^ Jo Jo 

By Point 3, Cf } (po) C -converges a.s. to c(po). Cg (po) tends to co(po) by assumption. 
(3 n (t) is a finite sum, so 

fin® - ^p.(po\p)Ci n) (p)^ n) (Po\p) 

p<p 

compactly. Finally, note that 

a n (s) =Ci n) ( Po ) (ao ^| n) («^» m ) + ^o ^ aC< n >(a,&,"0 

\ {a,b,m)dS (a,b,m)£S 

=Ci n \ Po )(a B^+b A^). 
By Lemma El and B[ n) converge compactly to 1 1— > 1/(1 + £), so 

lim a n (s) = ^-^-c s (p ) a.s. 

rw+oo 1 + t 

compactly. Since these are all compact convergences, we can pass to the limit in fl35l) . for 
all p G S. This readily shows that (c(p)) solves ( |T5|) . with initial conditions (co). 

□ 

Theorem 2. Assume the same hypothesis as in Proposition^ and assume as well 

(a 2 + 6 2 )co(a, 6, m) < +oo. 

(a,6,m)G5 

Lei T c 6e defined as in Definition^ Then {C^) t e[o.T c ) converges (in distribution) to the unique 
solution of Smoluchowski 's equation (pQ) . 

Remark 7. Obviously, convergence has to be understood with respect to Skorokhod's topology 
on [0, T) (which is the trace topology of Skorokhod's topology on [0, +oo)). In particular the 
sequence of the laws of (Cj; )te[o,T) is tight. 

Proof. Let Q n be the law of (C^)te[o,T c )- The sequence (Q n ) is tight. Let Q one of its limit 
points, and let c a process with law Q. By Proposition 0] above, c solves a.s. the system (|T5l) . 
with initial conditions (co). Now, let gt(x,y,z) the (a priori random) generating function of 
c. It is easy to see that g is well defined for (t,x,y,z) G [0,T C ) x (0, l) 2 x [0,1], and that 
g t (., .,z) is regular for every z G [0, 1]. Moreover, we see as in the proof of Proposition [3] that 
for every z G [0, 1], g t (., .,z) solves the PDE <HM with initial conditions (x, y) \— > g (x,y,z) = 
J2co(a,b,m)x a y b z m . Hence by Proposition [31 (q) solves Smoluchowski's equation [1] until T c . 
But by Theorem [U there is a unique solution to this equation on [0, T c ). Hence there is a unique 
limit point to (Q n ), so the sequence itself converges to the solution of Smoluchowski's equation 
on[0,T c ). □ 
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